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width  and  yields  models  equivalent  or  similar  to  the  traditional  Curtis-Godson 
approximation.  The  second  treats  approximations  to  the  spatial  derivatives  of 
these  functions.  From  the  standpoint  of  computing  line  or  band  radiance,  the 
spatial  derivatives  are  mere  fundamental  quantities  than  the  transfer  functions 
themselves;  consequently,  these  latter  "derivative  approximations"  arc  intrin- 
sically more  accurate  than  the  Curtis-Godson  type  approximations. 

All  of  the  models  are  formulated  to  give  the  spatial  derivative  of  the  mean 
equivalent  width  function  W(s)/6  in  the  form 

1 dW(s)  . . . . r.  , . . 

J ~ds~  = c(s)  p(s)  k(s)  y(s)  , 

[where  c(s),  p(s),  and  k(s)  are  the  concentration,  total  pressure,  and  absorption 
band  model  parameter,  respectively,  at  the  local  path  position  s.  The  func- 
tional form  of  the  derivative  function  y(s)  is  derived  for  each  of  the  models.  In 
general,  y(s)  is  a function  of  both  local  and  path-averaged  (subscript  e)  band 
model  parameters.  These  local  and  averaged_parameters  are  used  to  define  a ' 
dimcnsionlcjs  s_optical  depth  parameter  xe  = uke/0e  and  three  nonuniformity 
indexes  p y/ye,  r = B/Bt,»  and  q = 6e/*  (y  is  line  width,  6 is  the  mean  line 
spacing  band  model  parameter,  and  B 2rry/6).  The  explicit  manner  in  which 
these  parameters  enter  into  the  various  functional  forms  for  y(s)  is  derived  and 
discussed.  The  introduction  of  the  parameter  q is  an  important  aspect  of  the 
new  models  and  is  an  explicit  measure  of  nonisothe rmality  along  the  optical 
path.  In  addition  to  the  traditional  Curtis-Godson  definitiens_for  the  effective 
parameters  ke  and  Be,  new  definitions  for  the  path  averages  ye  and  ~E  are 
derived.  These  new  definitions  are  obtained  from  the  fundamental  properties  of 
the  assumed  inverse  line  strength  distribution. 

A summary  of  all  of  the  models  is  given  as  a table  of  relevant  equations  for  use 
in  practical  calculations. 
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I.  INTRODUCTION 


The  calculation  of  low- re  solution  radiance  and  transmittance  spectra 
of  molecular  gases  by  narrow  band  model  methods  has  widespread  applica- 
bility in  quantitative  spectroscopy.  Two  particularly  important  areas  of 
application  are  the  calculation  of  heat  transfer  in  the  earth's  atmosphere^  ^ 

(2 

and  the  calculation  of  infrared  radiance  signatures  of  hot  combustion  gases.  ‘ 
The  latter  spans  the  range  from  small  laboratory  flames  to  exhaust  plumes 
of  jet  aircrait  and  rocket-propelled  missiles. 

Most  practical  applications  of  band  model  methods  must  treat  nonuni- 
form optical  paths.  In  atmospheric  applications,  inhomogeneities  result 
from  the  variation  of  pressure  and  mixing  ratio  of  active  species  with  alti- 
tude. Temperature  variations  also  occur  in  the  atmosphere,  but  for  many 
applications,  these  variations  may  be  neglected  and  the  atmosphere  treated 
as  isothermal.  In  combustion  gas  applications,  however,  very  strong 
thermal  gradients  may  be  encountered.  This  is  particularly  ;rue  for  exhaust 
plume  problems  in  which  shock  structure  and  the  mixing  and  afterburning  of 
unburned  fuels  with  atmospheric  oxygen  are  considered.  An  important  prob- 
lem that  combines  aspects  of  atmospheric  and  combustion  gas  methods  is  the 
propagation  of  radiation  from  a hot  gas  source  through  a long  cool  absorbing 
atmosphere.  If  the  same  molecular  species  is  responsible  for  both  emission 
and  absorption,  the  high  degree  of  line  position  correlation  between  the  emis- 
sion and  absorption  spectra  may  preclude  the  decoupling  of  the  optical  path 
into  isolated  emitter  and  absorber  regions  and  multiplying  the  source  band 
radiance  by  the  absorber  band  transmittance  in  order  to  arrive  at  the  trans- 
mitted radiance  spectrum.  The  entire  optical  path  through  both  the  absorb- 
ing and  emitting  path  segments  must  be  treated  as  a single  nonuniform  optical 
path. 

Radiative  transfer  along  nonuniform  optical  paths  for  both  isolated  lines 
and  bands  of  lines  has  a long  and  varied  history  of  theoretical  study.  These 
studies  can  be  conveniently  grouped  into  three  categories  corresponding  to 


a priori  assumptions  concerning  the  type  of  optical  path  considered.  The 
first  category  includes  those  solutions  and  approximations  that  make  no  prior 
assumptions  concerning  the  variations  along  the  path  and  are  thus  applicable 
(but  not  necessarily  accurate)  to  any  general  nonuniform  path.  The  most 
significant  entries  into  this  category  include  the  Curtis-Godson  (CG)  ■'\pproxi- 
mation,  *3"  10)  the  Lindquist-Simmons  (LS)  approximation,*11  13^  tl  deriva- 
tive model  of  Lindquist  et  al.  , **^  the  intuitive  method  of  Weinreb  and 
Neundorffer,  *13)  and  the  exact  solution  for  an  isolated  Lorentz  line  of  Plass 
and  Fivel.*10)  The  formulation  discussed  and  applied  by  Simmons  and 
others*  l7~21)  can  also  be  considered  as  a general  path  approximation,  but 


requires  the  selection  of  parameters  or  functions  that  fit  the  model  to  spe- 
* (27' 

cific  types  of  nonuniform  paths.  The  work  of  Boughner  ‘ on  shock-heated 

gas  also  fits  into  this  restricted  general  path  category. 

The  second  category  includes  studies  that  consider  optical  paths  along 
which  the  pressure  is  specifically  assumed  to  vary  exponentially.  Applica- 
tion to  atmospheric  slant  paths  is  evident.  Significant  contributions  in  this 

(23)  • (24) 

category  include  the  works  of  Plass  and  Fivel,  Yamamoto  and  Aida, 

(25) 

and  Yamamoto  et  al. 

The  third  category  of  studies  are  most  applicable  to  the  problem  of 
flame  radiation  through  an  absorbing  atmosphere.  The  a priori  path  model- 
ine here  is  the  consideration  of  two  uniform  path  segments  in  series.  Treat- 
r*  ^ 2^  27 ) 

ment  of  line  radiation  for  this  case  has  been  made  by  Sakai  and  Stauffer 

[ with  comments  by  Plass*28)]  and  by  Goody.  *29)  Sakai  and  Stauffer  con- 
sidered the  general  problem  of  two  arbitrarily  overlapping  Lorentz  lines. 

For  zero  separation,  their  analysis  is  applicable  to  the  two-cell  problem. 
Treatment  of  band  radiation  for  two  path  segments  in  series  has  been  given 

by  Plass*30’3*)  and  Rodgers.*32) 

Despite  the  extensive  work  reflected  in  this  short  review,  no  band 
model  has  been  developed  that  adequately  handles  general  optical  paths  along 
which  high  degrees  of  inhomogeneity  and,  particularly,  nonisothermality 
occur.  The  most  often  used  of  these  approximations  is  the  CG  approxima- 
tion. Generally  speaking,  this  approximation  is  adequate  for  application  to 


» 
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atmospheric  problems  and  even  to  combustion  gas  problems  that  display 
moderately  strong  temperature  gradients.  1U’  For  very  strong 

thermal  gradients,  and  for  the  case  of  viewing  hot  gas  sources  through  inter- 
vening cool  absorption  paths,  the  CG  approximation  as  well  as  the  remaining 
general  path  formu1  tions  often  fail  to  yield  adequate  results. 

The  concern  of  the  present  paper  is  the  formulation  of  improved  band 
radiance  models  that  can  be  applied  efficiently  to  general  nonuniform  optical 
paths  and  that  account  for  large  nonuniformities  better  than  previous  models. 
Spectra  generated  with  these  models  and  experimental  spectra  are  compared 
in  a companion  paper/1  The  models  are  formulated  within  the  confines  of  the 
statistical  band  model  for  Lorentz  line  shape  and  inverse  line  strength  distri- 
bution. Three  works  are  of  particular  importance  in  establishing  the  ground- 
work for  these  formulations:  (1)  the  work  of  Godson^  ^ contains  the  original 
derivation  of  the  statistical  model  for  the  inverse  line  strength  distribution 
and  the  most  thorough  treatment  of  the  traditional  CG  approximation 
applied  to  the  statistical  model;  (2)  the  work  of  Lindquist  and  Simmons 
rationalizes  the  approach  of  approximating  the  derivatives  of  such  transfer 
functions  as  mean  transmittance  and  equivalent  width  rather  than  approxima- 
ting these  functions  directly;  and  (3)  the  recently  reported  work  of  Lindquist 
et  al/14^  provides  the  intuitive  direction  for  the  formulations. 

In  this  paper,  radiative  transfer  for  an  isolated  Lorentz  line  is  re- 
viewed in  Section  II,  and  the  statistical  band  model  for  uniform  optical  paths 
is  reviewed  in  Section  III.  Both  discu^ions  are  restricted  to  topics  required 
in  subsequent  developments,  but  are  presented  in  some  detail.  The  founda- 
tions established  in  these  two  sections  are  used  in  the  formulations  given  in 
Sections  IV  and  V.  In  Section  IV,  the  CG  approximation  as  given  by  Godson 
is  discussed,  and  an  alternative  formulation  of  this  approximation  is  con- 
sidered. This  alternative  formulation  establishes  interpretations  and 

' Stephen  J.  Young,  "Evaluation  of  Nonisothe rmal  Band  Models  for  I120, 

J.  Quant.  Spectr.  Radiative  Transfer  (to  be  published). 
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functional  relations  necessary  to  the  formulation  of  the  various  derivative 
approximations  made  in  Section  V.  The  latter  part  of  Section  IV  and  Sec- 
tion V comprise  the  significant  new  work  of  this  paper.  The  mo  e s are 
reviewed  and  summarized  in  Section  VI. 


II.  ISOLATED  LORENTZ  LINE 


A.  GENERAL  RELATIONS 

In  the  statistical  band  model  used  in  subsequent  formulations,  the 
average  transmittance  for  a spectral  interval  Av  is  expressed  in  terms  of 
the  integrated  absorptance  (equivalent  width)  of  the  individual  lines  m Av. 

In  this  section,  the  radiative  transfer  effects  of  an  isolated  Lorentz  line  are 
reviewed  to  the  extent  required  for  these  subsequent  formulations. 

Consider  an  optical  path  that  extends  from  the  geometric  position  s = 0 
^ s = S.  The  concentration  (mole  fraction)  of  active  gas  c(s),  total  gas  pre- 
sure  p(s),  and  temperature  T(s)  are  considered  as  general,  but  known,  func- 
tions of  s.  The  gas  is  assumed  to  be  in  a state  of  local  thermodynamic 
equilibrium.  The  spectral  absorption  coefficient  for  the  Lorentz  (dispersion) 

line  shape  is 


k(v,  s)  = 


S(s)  \(s) 


+ \"(s) 


where  S(s)  is  the  line  strength  and  \(s  s the  line  width  parameter.  The 
reference  of  spectral  position  is  the  line  center.  S(s)  is  assumed  to  be  a 
function  of  temperature  only,  whereas  y(s)  may  be  a function  of  all  three  of 
the  path  variables.  The  variations  of  S(s)  and  \(s)  with  c(s),  p(s),  and  T(s) 
are  presumed  known. 

The  two  fundamental  quantities  of  line  radiation  that  are  important  here 
are  the  line  radiance  and  the  line  ab  so-ptance.  The  former  is  define  d by 


N 


(1) 


where  W(t)  is  the  integrated  absorptance  of  the  line  (equivalent  width)  for 
the  peth  segment  between  s = 0 and  the  general  position  s and  is  defined  by 


W{s) 


c(s/)p(s/)S(s/)y(s/) 

2 , 2.  / . 
v + y (s  ) 


(2) 


The  line  absorptance  is  the  equivalent  width  evaluated  for  the  whole  path  W(S). 
N*(s)  is  the  Planck  radiation  function  evaluated  at  the  true  spectral  position 
of  the  line  center  and  at  the  temperature  T(s).  The  form  of  equation  (1)  is 
based  on  the  assumption  that  N:':(s)  changes  negligibly  over  the  spectral  width 
of  the  line. 

B.  EXACT  SOLUTIONS  FOR  W(s) 

The  prediction  of  either  N or  W(S)  for  the  line  reduces  formally  to  the 

evaluation  of  W(s)  for  general  s.  For  a uniform  optical  path  along  which  c(s), 

(34) 

p(s),  T(s),  and  hence  S(s)  and  •y(s)  are  constants,  the  solution'  to  equa- 
tion (2)  is 


W(s)  = ZvyL  [x(  s)]  , (3) 

where  x(s)  = u(s)S/27i-y  is  a dimensionless  optical  depth  parameter,  and 
u(s)  = cps  is  the  optical  depth  between  s - 0 and  s.  L(x)  is  the  Ladenburg- 
Bei  he  function 


L(x)  = xe"X[IQ(x)  + I ^ ( x )]  , (4) 

where  lf.(x)  and  I^(x)  are  modified  Bessel  functions.  L(x)  is  shown  in  Fig.  1. 
The  asymptotic  limits  of  L(x)  for  weak  (x  - 0)  and  strong  (x  — ®)  absorption 
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Figure  1.  Curve -of-G  rowth  Functions.  — Landenburg-Reiche 

functions  L(x); inverse  distribution  function 

I(4x)/4; exponential- tailed  inverse  distribution 

function  I ( ttx)  / it. 


are  L(x)  — x and  L(x)  — yj 2x/ tt,  respectively.  The  corresponding  limits  for 
W(s)  are 


W(s) 


u(s)S  = 2t ryx 
2%/u{s)Sy  - 2^yJ~ 


weak 


strong 


(5) 


An  exact  solution  for  W(s)  for  an  arbitrary  opticaJ  path  has  been  given 
by  Plass  and  Fivel.  The  solution  is  a series  expansion  in  terms  of 

moments  of  ■y(s).  The  complexity  of  the  solution  and  the  fact  that  only  the 
first  three  expansion  coefficients  have  been  evaluated  preclude  its  use  in 
most  practical  applications.  This  solution  is  not  considered  further  here. 

Jn  the  limits  of  weak  and  strong  absorption,  general  solutions  of  equa- 
tion (2)  can  be  obtained  in  simple  form.  Familiar  arguments^  used  to  ob- 
tain the  weak  and  strong  limits  for  a uniform  path  are  readily  applied  to  a 
nonuniform  path  and  yield 


If  c(s')p(s')S(s')  ds' 

2 y J c(  s ' )p(  s ' )S(  s ' )y(s  ' ) ds' 


weak 


strong 


At  this  point,  it  is  convenient  to  introduce  the  concept  of  path-averaged 
line  parameters.  The  path-averaged  line  strength  S^(s)  and  line  width  -VC1S) 
are  defined  by 


S (s) 
e 


1 

u(  s) 


c( s ' )p( s ' )S( s') 


d s ' 


(6) 


-iZ- 


and 


Ve(s)  = u(s)SoTO  f0  cU'lpU'JSU'WU'td.' 


(7) 


where  u(s)  is  the  optical  depth 


u(s)  = f c(s'  )p(s ' ) ds ' 

Ja 


(8) 


A convenient  dimensionless  optical  depth  parameter  is 


u(s)Se(s) 

se<s)r  2ne('sT 


(9) 


In  terms  of  these  path-averaged  parameters,  the  limiting  solutions  for  W(s) 
may  be  written  as 


W(s) 


u(s)Se(s)  = 27T^e(s)xe(s) 



2s/u(s)Se(s)ye(s)  = 2iye(s)yJ 


weak 


strong 


(10) 


Comparison  of  equations  (9)  and  (10)  indicates  that  W(s)  for  a general 

path  approaches  the  same  limiting  forms  for  weak  and  strong  absorption  as 

does  W(s)  for  a uniform  path.  The  difference  is  that  the  path  averages  Se(s) 

and  v (s)  replace  the  constants  S and  y,  respectively,  and  the  general  result 
' e 

for  u( s ) given  by  equation  (8)  replaces  the  quantity  cps. 


Considerations  presented  later  require  the  spatial  derivative  of  these 
exact  limiting  solutions.  Straightforward  differentiations  yield 


dW(  s) 
ds 


c(s)p(s)S(s) 


< 


c( s)p( s)S( s)y(s) 
yU(s)Se(s)ye(s) 


c( s )p( s )S( s)  \(s) 


weak 


strong 


(ID 


C.  CURTIS -GODSON  APPROXIMATION 

(3)  (4,  5) 

The  CG  approximation  was  introduced  by  Curtis  and  Godscn  ' in 
the  context  of  isothermal  band  models  used  for  computing  heating  rates  in 
the  atmosphere.  An  extensile  literature  exists  on  the  application  and  inter- 
pretation of  this  approximation.  The  approximation  can  be  formulated  from 
several  points  of  view'.  The  conceptually  simplest,  and  the  one  most  in 
keeping  with  Curtis'  original  discussion,  is  that  the  path- integrated  line- 
shape  function  [the  exponential  argument  of  W(s)  in  equation  (2)]  is  replaced 
by  a Lorentz  line- shape  function  that  depends  on  effective  line  parameters 
S(s)  andy(s).  That  is,  W(s)  is  approximated  by 


W(s) 


u(s) 


S(s) 


TT 


2 

V 


5ii) | 

+ y2(s)Jj 


dv 


An  immediate  Ladenburg-Reiche  integration  can  be  made  to  give 


W(s)  = 2tt^L[x(  s)] 

Aith 


x(s) 


u(s)S( s] 
2iry(  s ) 
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The  effective  parameters  S(s>  and  y(s)  are  determined  by  forcing  the 
approximate  solution  for  W(s)  to  agree  with  the  exact  solutions  [equation  (10)] 
for  weak  and  strong  absorption.  The  result  is  that  the  effective  parameters 
arc  best  defined  as  the  path-averaged  parameters  of  equations  (o)  and  (7), 
i.  e.  , S(s)  = Se(s)  and  y(s)  = V e(»>.  The  CG  approximation  for  an  isolated 
Lorencz  line  and  general  nonuniform  optical  path  is,  then, 


W(s)  = 2iry  ( s)L[x  ( s)l 


(12) 


The  functional  form  is  the  same  as  for  a uniform  path,  but  the  constant 
parameters  S and  y are  replaced  by  the  path  averages  Se(s)  and  \£(s), 

respectively. 

The  calculation  of  line  radiance  according  to  equation  (1)  requires  the 
derivative  dW(s)/ds.  Differentiation  of  equation  (12)  yields 

dWU)  _ c(s)p(s)S(s)y[x  (s),  p(s)]  , ^ 

QS  e 


where 


y(x,  p)  = (2 


r dL(x) 

P)_dT" 


+ (P  - l) 


L(x) 


(14) 


with  dL(x)/dx  - e"xl0(x).  The  parameter  p(s)  is  the  ratio  o 
of  line  width  to  the  path-  averaged  line  width. 


f the  local  value 


i(s)  (15) 

\ev> 


and  is  the  index  that  measures  the  degree 


of  local  path  nonuniformity. 
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The  asymptotic  limits  for  the  derivative  function  y(x,  p)  given  by 
equation  ( 14)  are  readily  verified  to  be  y(x,  p)  -»  1 as  x - 0,  and  y(x,  p)  -» 
p /y/Zttx  as  x - 00 . Comparison  of  equation  (13)  with  the  exact  solutions  of 
equation  (11)  reveals  that  these  limits  for  y(x,  p)  are  the  correct  limits. 
Although  y(x,  p)  approaches  the  correct  limits  for  weak  and  strong  absorp- 
tion, spurious  behavior  may  occur  for  intermediate  values  of  x if  the  non- 
uniformity index  p(s)  is  large.  A large  value  of  p(s)  would  occur,  for 
example,  if,  along  an  isothermal  path,  there  were  a sudden  pressure  in- 
crease such  tl  it  y(s)  » y^(s).  The  derivative  functii  n is  shown  in  Fig.  2 
as  a function  of  x for  a few  values  of  p.  It  can  be  shown  that  if  p > 2,  there  is 
a region  of  x values  for  which  y(x,p)  > 1.  From  equation  (11),  it  is  evident 
that  this  effect  predicts  that  W(s)  grows  at  a race  faster  than  does  for 
weak  absorption.  This  behavior  is  unacceptable  and  points  out  a funda- 
mental flaw  in  the  CG  approximation. 

D.  LINDQUIST- SIMMONS  APPROXIMATION 

(11' 

The  approximation  of  Lindquist  and  Simmons'  ' corrects  the  flaw  in 
the  CG  approximation  noted  above.  This  approximation  proceeds  from  the 
observation  that  the  line  radiance  given  by  equation  (1)  is  formulated  in 
terms  of  the  derivative  dW(s)/ds  rather  than  W(s)  directly.  The  LS  approxi- 
mation thus  considers  a direct  approximation  to  dW(s)/ds  lather  than  to 
W(s).  The  spatial  derivative  of  W(s)  from  equation  (2)  is 


dW(s) 

ds 


c(s)p(s)S(s) 


\2(s) 


exp 


c(s/)p(s/)S(s/)\(s/)  , 

2,2,/.  d& 
v + \ (s  ) 


d v 


As  ii  the  CG  approximation,  the  path-integrated  line- shape  function  (the 
exponential  argument.)  is  replaced  by  a Lorentz  line- shape  function  expressed 
in  terms  of  effective  line  parameters  S(s)  and  \(s).  The  derivative  may  then 
be  expressed  as 


= c(s)p(s)S(s)y[x(s),'p(s)] 


(16) 
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■“J 


whe  re 


and 


x(s) 


u( s)S( s) 
2tty(s) 


\(s) 


y(x 


n 


( 2x  ) dz 
exp  j < 


1 + z2p2J  1 + z2 


(17) 


A simple  analysis  of  equation  (17)  yields  the  asymptotic  limits  of  y(x,  p)  as 

y ( x:,  p)  - 1 as  x - 0,  and  y(x,  p)  - p / J It rx  as  x — Thus,  the  weak  absorption 

limit  of  dW(s)/ds  given  by  equation  (16)  reduces  to  the  correct  limit  given  by 

equation  (11)  regardless  of  the  choice  of  S(s)  and  y(s).  Agreement  between 

***•» 

equations  (16)  and  (11)  for  strong  absorption  requires  S(s)\(s)  = S (s)\  (s), 

C G 

which  can  be  satisfied  with  the  intuitive  choices  S(s)  = S^fs)  and  y(s)  = y^(s). 
The  final  result  for  the  LS  approximation  is 


~ c(s)p(s)S(s)y[xc(s),  p(s)] 


(18) 


An  alternative  approach  to  the  CG  approximation  is  to  replace  \(s/)  in  only 
the  denominator  of  the  exponential  argument  of  equation  (2)  by  an  effective 
w idth  parameter  y(  s).  A Ladenburg-Reiche  integration  and  a forced  fit 
between  the  approximate  and  exact  solution  for  W(s)  in  the  weak  absorption 
limit  yields  \(s)  = Ye(s)  and,  finally,  the  result  of  equation  (12).  If  this  same 
ipproach  is  used  in  the  LS  approximation,  the  correspondence  \(  s)  = \ (s) 

G 

is  uniquely  determined  by  the  strong  absorption  fit  of  the  approximate  to 
exact  solution  for  dW(s)/ds. 
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The  derivative  function  of  equation  (17)  is  not  expressible  in  terms  of 

elementary  functions.  A tabulation  and  approximate  solutions  for  extreme 

(12) 

values  of  the  arguments  x and  p have  been  given  by  Young.  The  function 

is  shown  in  Fig.  2 for  comparison  with  the  CG  approximation.  The  function 
is  well-behaved  and  never  exceeds  unity.  Since  the  LS  approximation  treats 
dW(s)/ds,  no  solution  for  W(s)  is  given  other  than  the  formal  solution 


W(s) 


dW(s') 

ds 


ds' 


which  can,  in  application,  be  evaluated  numerically. 
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III.  STATISTICAL  BAND  MODEL 


A.  GENERAL  RELATIONS 

In  this  section,  the  radiative  transfer  effects  are  reviewed  for  an  array 
(band)  of  many  lines,  all  with  centers  confined  to  a spectral  interval  Av. 

The  average  radiance  in  Av  due  to  these  lines  and  for  a general  nonuniform 
optical  path  is 


N* ( s ) ds  . (19) 

d s 

This  form  of  the  band  radiance  equation  requires  Av  to  be  small  enough  that 
the  Planck  function  N*(s)  can  be  considered  as  sensibly  constant  across  Av. 
t(s)  is  the  average  spectral  transmittance  for  Av  and  the  path  segment 
between  s = 0 and  general  position  s.  The  exact  expression  for  t(s)  is 

T(S)  = A ~f  T(v»s)dv  » (2°) 

Av 

where  t(v,  s)  is  the  spectral  transmittance  within  Av.  If  the  number  of  lines 
N in  Av  is  large,  t(v,  s)  is,  in  general,  a rapidly  fluctuating  function  of  v, 
and  the  evaluation  of  t(s)  by  equation  (20)  is  quite  tedious.  In  addition  evalua- 
tion of  t(s)  by  this  method  requires  detailed  knowledge  of  the  spectral  posi- 
tion, strength,  and  width  of  each  of  the  N lines  in  Av.  In  order  to  expedite 
the  evaluation  of  t(s),  band  model  techniques  are  used.  The  objective  of 


All  of  these  lines  are  assumed  to  arise  from  a single  molecular  gas  species. 
Treatment  of  more  than  one  species  is  trivially  handled  by  invoking  the 
multiplicative  property  of  band  transmittance  for  nonidentical  species. 


Preceding  page  blank 
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band  model  methods  is  to  approximate  t(s)  by  an  expression  that  depends  on 
only  a feu  parameters.  These  band  model  parameters  reflect  average 
pioperties  of  the  lines  in  Av. 

The  statistical  band  model  approximation  of  t(s)  is  developed  from  two 
fundamental  assumptions.  One  is  that  the  placement  of  lines  in  Av  can  be 
modeled  by  a random  arrangement  of  equal  width  lines.  1 he  result  of  this 
assumption  is  that^ 


. T W(s)l 

t(s)  = exp  - — — 


where  W(s)  is  the  average  value  of  W(s)  for  the  lines  in  Av, 


(21) 


N 

W(s)  = ^ £ Wti.s)  * (22) 

i=  1 

and  1 is  a mean  line  spacing  parameter  6 = Av/N.  This  step  in  the  formula- 
tion eases  the  difficulty  in  evaluating  t(s)  to  the  extent  that  the  individual 
positions  and  widths  of  the  lines  in  Av  are  no  longer  required.  The  strengths 
are  still  required  for  computation  of  W(i,  s). 


This  expression  for  t(  s)  is  the  limiting  form  for  N -*  <®  and  Av  -»  °°.  The 
latter  ensures  that  the  absorption  caused  by  each  line  in  Av  is  confined  to 
A\ . An  additional  note  of  importance  is  that  r(s)  derived  in  the  statistical 
model  is  not  a mean  spectral  transmittance  in  the  sense  of  equation  (20), 
but  rather  a spectral  transmittance  (for  some  position  v in  Av,  generally 
the  center  of  A^)  that  has  been  averaged  over  all  possible  line  arrange- 
ments. It  can  be  established,  however,  for  the  condition  Av  that 
equation  (21)  results  for  any  choice  of  reference  position  v in  Av.  Thus, 
equation  (21)  is  also  an  adequate  representation  for  the  average  spectral 
t ransmittance. 
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The  other  modeling  assumption  is  that  the  strengths  of  the  lines  in  Av 
can  be  described  by  a continuous  distribution  function  P(S).  Then,  rather 
than  writing  the  average  value  of  W(s)  as  a summation,  as  in  equation  (22), 
the  integral  expression 


W(s) 


r 


P(S)W(S,  s)  dS 


(23) 


is  used.  The  use  of  this  definition  requires  the  se^ction  of  a plausible 
distribution  function  P(S)  and  the  expression  of  W(s)  for  an  isolated  line  in 
terms  of  a line  strength  parameter  S.  For  an  isothermal  path  (where  S is  a 
constant),  this  latter  requirement  is  trivial.  For  a noni sothermal  optical 
path,  the  expression  of  W(s)  [and  its  ch  rivative  dW(s)/ds]  in  terms  of  an 
acceptable  strength  parameter  is  the  basic  problem  considered  in  Sections  IV 
through  VI.  For  the  remainder  of  this  section,  a uniform  optical  path  is 
assumed. 

B.  INVERSE  LINE  STRE  GTH  DISTRIBUTION 

P(S)dS  is  the  probability  that  an  arbitrarily  selected  line  in  Av  will 

have  strength  S in  the  interval  dS  about  S.  Several  distribution  functions  have 

been  considered.  Use  of  the  dirac  distribution  P(Sj  = 6(S  - S)  is  equivalent 

to  the  assumption  that  all  lines  in  Av  have  a strength  equal  to  the  average 

line  strength  S.  The  use  of  an  exponential  distribution  P(S)  = S exp(-S/S) 

(35  36) 

yields  the  Mayer-Goody  form  of  the  statistical  model.  ’ ' The  distribu- 

tion functions  used  here  are  the  truncated  inverse  line  strength  distribution 

introduced  by  Godson'  ’ ’ and  the  exponential-tailed  inverse  distribution 

( 38 ) 

introduced  by  Malkmus.  T'ae  use  of  these  two  distributions  yields  nearly 

identical  results  for  man,  applications.  In  the  present  work,  the  distribu- 
tions are  used  interchangeably  in  the  following  manner:  The  formulation  of 
models  is  made  with  the  inverse  distribution;  the  functions  required  to  apply 
the  models  are  derived  with  the  exponential-tailed  inverse  distribution. 
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The  form  of  the  inverse  distribution  used  here  is 


o o < s < sm/r 


P(S)  = \ 


1 


S In  R 
0 


SM/R  " S £ SM 


(24) 


S > S 


M 


where  SV1  and  SM/R  are,  respectively,  the  maximum  and  minimum  line 
strengths  allowed  in  the  distribution  and  are  introduced  so  that  P(S)  can  be 
normalized.  SM  and  R are  two  independent  parameters  of  the  distribution. 
R is  generally  assumed  to  be  large.  The  mean  line  strength  for  the  distri- 
bution is 


S = 


R - 1 
R In  R 


* 


with  S - SM/ln  R for  large  R. 

For  application  to  nonisothermal  paths,  it  is  necessary  to  know  in 
more  explicit  form  how  temperature  enters  into  the  distribution.  In  order 
to  accomplish  this,  a functional  relationship  between  line  strength  and 
temperature  must  be  assumed.  The  form  assumed  here  is 


S(i,  T)  - SM(i,T) 


- 0 ( i ) / T 
e 


(25) 


where  i is  a line  index,  and  0(i)  is  the  energy  (expressed  in  temperature 
units)  of  the  lower  level  of  the  quantum  transition  causing  the  line.  This 
form  for  S(i,  T)  follows  directly  from  the  assumption  that  the  strength  in 
absorption  is  the  product  of  the  number  of  molecules  in  the  lower  level  and 
an  intrinsic  oscillator  strength  CT(i).  A dominant  factor  that  influences  the 
magnitude  of  o(i)  is  the  energy  change  in  the  transition  causing  the  line. 
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lines  confined  to  a small  spectral  interval  Av,  this  energy  change  is  nearly 
constant.  To  a first  approximation,  then,  a(i)  can  be  considered  independent 
of  i.  The  preexponential  factor  SM(i,  s)  then,  is  also  independent  of  i.  This 
simplification  is  used  and  yields  the  final  assumed  form  for  the  temperature 
variation  as 


a,-  -0(i)/T 

S(i,  T)  = ) e 


(26) 


This  form  defines  a unique  relation  (for  fixed  T)  between  the  strength  of  a 
line  and  the  energy  of  the  lower  level  of  the  transition  causing  the  line.  A 
distribution  of  line  strengths  thus  implies  a distribution  of  energy  levels. 

The  use  of  the  inverse  strength  distribution  of  equation  (24),  the  relationship 
between  S and  0 given  by  equation  (26),  and  the  transformation  relation 
P(S)d(S)  = P(0  )d(0  ) gives  the  energy  level  distribution  as 


P(0)  = { 


0 

1 


T In  R 

0 


0 < 0 

O<0  5 T In  R 
0 > T In  R 


Because  the  strongest  lines  in  Av  occur  for  absorption  from  (or  emission  to) 
the  ground-level  0=0,  S^(T)  in  equation  (26)  can  be  identified  as  the  distri- 
bution parameters  in  equation  (24).  The  distribution  P(0  ) is  constant  in 
the  variable  0 up  to  a maximum  value  0m  = T In  R.  Since  the  distribution  of 
energy  levels  within  a molecule  is  an  intrinsic  property  of  the  molecule,  we 

must  insist  that  0 be  independent  of  temperature.  This  consideration  thus 
ni 

yields  an  explicit  form  for  the  temperature  dependence  of  the  distribution 
parameter  R as 


R(T ) 


0 

e 


m 


/T 
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The  assumption  that  R is 
the  inverse  line  strength 


large  implies  that  0jn  is 
distribution  used  here  is 


large.  The  final  form  of 
(with  R and  Gm  " 00 ) 


P(S,  T) 


T 

c 5 


m 


0 


0 - S - 

S > VT) 


(27) 


where  the  explicit  dependence  on  temperature  is  included. 
and  3 are  considered  as  the  tVk°  indePendent  parameters 

C,  STATISTICAL  MODEL  FOR  A UNIFORM  PATH 


and  where  SM(T) 
of  the  distribution. 


For  an  isolated  Lorentr  line  and  a uniform  optical  path,  W is  Riven  as 
a function  of  S throueh  the  Ladenburg-Reiche  solution  [equation  (3)]  W 
2rvL(x),  a - Su/2n.  From  equations  (23)  and  (27),  the  mean  equivalent 
width  function  is  approximated  as  (with  the  transformation  7,  = Su/2^) 


w 

6 


dT] 


where 


x 


M 


S..U 

M 

2rv 


‘ When  the  distribution  is  written 
if  the  inverse  distrioution  is  as 
holds  for  all  temperatures. 


in  this  form,  it  can  be  readily  verified  that 
sumed  to  hold  for  any  temperature,  then  it 


a 
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The  integral  over  Tj  is  easily  evaluated  by  using  equation  (4)  for  L(  ,)  and 
standard  tables  of  integrals  for  Bessel  functions.  The  result  for  the  integral 

is 

I( x ) = f Mp  dT]  = 2xe"X[lo(x)  + I^x)]  -f  e XIQix)  - 1 . (28) 

*'0 

Thus,  the  solution  for  W/6  is 


) 


Introduction  of  the  parameters 


- SMT  _ S 

k " W"  c 


m 


1 _ 4T 

T 60 
6 m 


3 = 


2ttv 


and 


x 


ku 

= 3 


permits  the  result  to  be  simplified  to 


W 

6 


= | I(4y  ) 
4 


(29) 


(30) 


(31) 


(32) 


(33) 
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Equations  (28)  and  (33)  constitute  the  result  first  obtained  by  Godson.15' 37  > 
The  curve-of-growth  function  I(4x)/4  is  shown  in  Fig.  1.  Thewea  an 
strong  absorption  limits  are  the  same  as  for  the  Ladenburg-Reiche  function, 
but  the  range  of  x over  which  the  transition  occurs  is  broader.  This  phe- 
nomenon  occurs  because  the  st.ong  aosorp.ion  limit  for  a band  of  lines  >s  no, 
obtained  until  a sizable  fraction  of  the  weakest  lines  allowed  by  the  line 
strength  distribution  display_stvong  absorption  behavior.  The  weak  and 
strong  absorption  limits  of  W/6  are 


W 

6 


uk 


w 


eak 


strong 


D>  BAND  MODEL  PARAMETERS 

Equation,  (29)  through  (31)  relate  the  two  parameters  k and  0 to  the 
arrangement  parameter  6 and  the  strength  distribution  parameters  bM 
-ind  9 . These  utter  three  parameters  have  meaning  only  within  the  con- 

fines of  the  modeling  assumptions  that  lead  to  the  results  of  equation  (33) 

ThL.  practical  interpretation  of  band  modeling  procedures  is  that  the  resu  ts 
(equations  (28)  and  (33)  in  this  case]  constitute  a functional  form  for  the 
variation  [by  equation  (2,)]  with  optical  depth  u that  depends  on  two 

band  model  parameters  k and  3 that  are  to  be  determined  empirically  by 

fitting  the  functional  form  to  experimental  data. 

Two  general  procedures  are  used  to  determine  k and  0.  In  one,  ex- 
perimental absorption  measurements  are  made  on  a uniform  gas  sample  at 
fixed  temperature  V with  a spectrometer  centered  at  the  spectral  position 
of  interest  and  adjusted  to  a spectral  resolution  A..  The  optical  thickness  u 
„f  the  sample  is  varied  over  a wide  range  such  that  both  weak  and  strong 

pi  ion  .lata  are  obtained.  By  fining  the  weak  absorption  portion  of  the 

experimental  curve  7(u>  to  the  model,  k can  be  determined.  By  fitting  the 
strong  absorption  data,  the  product  k0,  and  hence  0,  can  be  obtained.  n 
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this  manner,  the  parameters  k and  3 can  be  obtained  for  all  spectral 
intervals  Av  across  an  entire  molecular  absorption  band  and  for  a wide  range 
of  temperatures. 

'.Che  other  approach  to  the  determination  of  k and  3 can  be  used  if  the 
line  optical  parameters  S(i)  and  ■y(i)  are  known  for  all  of  the  lines  in  each 
interval  Av.  The  temperature  dependence  of  S(i)  and  yU)  must  also  be  known, 
From  equations  (22)  and  (3),  W/6  for  a uniforn  optical  path  may  be  written 
in  terms  of  the  line  parameters  S(i)  and  y(i)  as 


W 

6 


(3  0 


Equation  (33)  represents  this  same  quantity  expressed  in  terms  of  the  band 
model  parameters  k and  3.  By  forcing  equations  (33)  and  (34)  to  be  identical 
for  weak  and  strong  absorption,  the  following  definitions  can  be  obtained: 


i=  1 


N 


0 


N 


t2 


£ v/S(i)^(i) 

i=  1 J 


(35) 


(36) 


In  both  this  and  the  experimental  absorption  approach,  k and  3 are  deter- 
mined  by  fitting  the  model  in  the  limits  T = 0 and  t = ».  Godson  ' has 
indicated  that,  for  some  applications,  fits  at  other  transmittance  values 
may  be  desirable.  In  his  work  on  the  atmospheric  transmittance  of  H^O, 
for  example,  he  used  fits  at  t - 0.  45  and  0.  95  to  obtain  a best  overall 
agreement  between  model  and  data. 
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Nolo  that,  by  including  -y(i)  in  the  definition  for  & as  a parameter  that  depends 
on  the  line  index  i,  the  model  approximation  that  all  -y  are  constant  has  been 
lessened. 

Although  only  the  two  parameters  k and  3 are  required  for  application 
of  the  isothermal  band  model,  it  is  convenient  (and  necessary  for  some  ex- 
tensions of  the  band  model  formulation  to  nonisothermal  paths)  to  consider  3 
as  a ratio  of  two  parameters  y and  > by  a generalization  of  equation  (31), 


3 = 12Y 

6 


(37) 


where  y is  a measure  of  the  average  line  width  for  the  lines  in  Av,  and  6 is  a 
measure  of  the  mean  line  spa  .»ig  in  Av.  The  use  of  three  parameters  (k,  6, 
and  y)  in  a two-parameter  model  required  auxiliary  information  on  one  of  the 
parameters.  This  information  is  generally  introduced  by  performing  inde- 
pendent experimental  measurements  of  y or  by  making  theoretical  and  em- 
pirical assumptions  on  y.  When  band  model  parameters  are  determined 
from  individual  line  data  [equations  (35)  and  (36)],  y is  conveniently  defined 
by 


V = 


N 


N 


Z VU) 
ml 


An  advantage  of  introducing  the  two  parameters  y and  7 in  place  of  6 is  that 
the  variations  of  3 with  the  composition  of  the  gas,  i.  e.  , the  ratio  of  active 
gas  concentration  to  various  foreign  gas  concent  rations,  can  all  be  placed  in 
which  leaves  6 as  a function  of  temperature  (and  spectral  position)  only. 
When  the  fundamental  band  model  parameters  "k,  7,  and  y have  been  deter- 
mined as  a function  of  spectral  position  and  temperature,  it  is  possible  to 
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define  the  modeling  parameters  and  the  product  66^.  From  equations  (29) 
and  (30),  the  following  expressions  can  be  obtained: 


Sm(T)  = 4k(T)6(T) 


(38) 


66  = 4T6(T) 

m 


(39) 


These  two  relations  are  particularly  important  in  the  formulation  of  ncn- 
isothermal  band  models  in  Sections  IV  and  V.  An  interesting  side  aspect 
of  equation  (39)  is  that  the  line  density  parameter  Ds  1/6  i. s predicted  to 
display  a linear  temperature  dependence.  This  result  offers  an  opportunity 
to  experimentally  determine  the  adequacy  of  the  entire  modeling  procedure 
by  examining  the  constancy  of  the  product  T6  with  temperature. 

E.  EXPONENTIAL-TAILED  INVERSE  LINE  STRENGTH 
DISTRIBUTION 

Same  comments  on  the  use  of  the  exponential-tailed  inverse  line 

strength  distribution  rather  than  the  inverse  distribution  are  appropriate 

(38) 

here.  This  distribution  function  is' 


P(S)  = 


_1 

S In  R 


r -s/sm 


-RS/S./I 

M 

e 


where  S is  allowed  to  assume  all  values  Os  S £ ®.  is  no  longer  inter- 

preted as  a maximum  allowed  line  strength,  but  simply  as  a parameter  of 
the  distribution.  The  distribution  results  by  use  of  the  more  general  result 
of  equation  (25)  over  (26)  for  line  strength  and  the  assumption  that  sM(i,  T) 

is  distributed  about  S.  ,(T)  according  to  an  exponential  strength  distribution. 

iVl 

In  a sense,  it  is  more  general  than  the  purely  inverse  line  strength 
distribution. 
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For  a uniform  optical  path,  the  formulation  of  a band  model  with  the 

exponential-tailed  distribution  parallels  the  formulation  with  the  inverse 

distribution,  except  that  the  curve-of-growth  function  I(4x)/4  is  replaced 

1/2 

by  I ( ttx ) / tt,  where  Mx)  = (1  + 2x)  - 1,  and  the  factor  4 in  equations  (30), 

(33),  (38),  and  (39)  is  replaced  by  tt. 


IV.  CURT  IS -GODSON  APPROXIMATIONS 


In  this  and  Section  V,  the  fundamentals  reviewed  in  Sections  II  and  III 
are  used  to  formulate  band  models  for  general  inhomogeneous,  nonisothermal 
optical  paths.  Here,  two  approaches  to  the  familiar  CG  approximation  for 
bands  of  lines  are  considered:  that  originally  given  by  Godson^  and  one 
that  is  based  on  the  same  assumptions  used  in  Godson's  approach,  but  that 
expands  on  the  nature  of  the  assumptions.  The  second  appioach  also  yeilds 
expressions  that  are  essential  to  the  formulations  in  Section  V. 

A.  GODSON'S  FORMULATION 

In  Godson's  approach,  we  assume  that  the  functional  form  of  the  curve 
of  growth  for  a nonuniform  path  is  the  same  as  that  for  a uniform  path.  The 
expression  for  the  mean  equivalent  width  is  written,  from  equation  (33),  as 


W(s) 

6 


3e(s) 

4 


I[4xe(s)] 


u(s)k  (s) 
xe(s)  - Be(s) 


(40) 


where  k (s)  and  8 (s)  are  effective  band  model  parameters  to  be  determined, 
e e 

and  the  optical  depth  u(s)  is  given  by  equation  (8).  The  effective  parameters 
k (s)  and  0 (s)  are  determined  by  a procedure  similar  to  that  used  to  deter- 
mine k and  0 in  terms  of  line  parameters  for  a uniform  path.  That  is,  the 
weak  and  strong  absorption  limits  of  equation  (40)  are  forced  to  agree  (as 
well  as  possible)  with  the  w'eak  and  strong  results  obtained  from  W(s)/6  given 
by  equation  (22)  and  the  exact  limiting  solutions  to  W(i,  s)  given  by  equation 
(10).  The  resulting  definitions  for  k^s)  and  0g(s)  parallel  the  definitions  of 


equation  (35)  for  k and  equation  (36)  for  3,  except  that  they  are  written  in 
terms  of  the  path-averaged  parameters  S^(i,  s)  and  \e(i,  s).  They  are 


N 

ke(s)  = SS  E 
i=  1 


S£(i,  s) 


(41) 


and 


3 „( s ) = 

e 


2tt 

1 

ke(s) 

Av 

N 


i=  1 


12 


(42) 


These  definitions  are  not  particularly  useful  because  they  not  only  require  a 
detailed  knowledge  of  individual  line  parameters,  but  also  require  that  the 
path  averages  of  the  strength  and  width  of  each  line  be  recomputed  for  each 
nonuniform  path  considered.  These  relations  can  be  used,  however,  to 
formulate  definitions  that  are  practical.  The  procedure  is  most  easily 
shown  for  k^ts).  By  interchanging  the  order  of  line  summation  and  path  inte- 
gration implicit  in  the  above  definition  for  k^(  s ),  we  obtain 


k (s) 
e 


1 

u(s) 


r s 

N 

/ c( s ' )p(s ' ) 

i;Ls(i's/) 

'o 

i=  1 

d s ' 


From  equation  (35),  the  quantity  in  brackets  can  be  identified  as  the  parame- 
ter k that  corresponds  to  an  isothermal  path  at  temperature  Tfs7).  Conse- 
quently, 


k (s) 
e 


c(s')p(s')k(s') 


ds ' 


(43) 


Sp(i,  s)  and  VpU>  s)  are  defined  by  trivial  generalizations  of  equations  (6)  and 
(7),  respectively,  that  include  the  line  index  i as  an  explicit  argument  to  S 
and  \. 
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where  k(s')  denote s k[T( s')].  The  effective  absorption  band  model  parameter 
k (s)  is  thus  the  average  value  of  k over  the  nonisothermal  path.  Ifjc  is 
known  for  all  temperatures  in  the  range  that  occurs  along  the  path,  k^s)  can 

be  computed  for  all  path  positions  s. 

In  order  to  obtain  a corresponding  result  for  3e(s),  an  additional  as- 
sumption must  be  made  in  order  to  permit  the  interchange  of  the  order  of 
integration  and  summation  in  equation  (42).  With  the  definitions  of  equations 
(6)  and  (7),  0 (s)  may  be  written  as 


The  additional  assumption  made  by  Godson  is  that  S(i,  s ) and  y(  . s ) may 
be  written  in  the  function*.!  product  forms 


S(i,  s')  = a^Oa^s') 


(44) 


and 


y(i,  s')  = gj(i)g2(s'  ) . (45) 

The  use  of  these  relations  and  an  interchange  of  the  order  of  summation  and 
integration  yields 

0 (s)  ^ fS  c(s')p(s') 

e u(s)k  (s)*q) 

The  quantity  under  the  square  root  sign  is  identified,  from  equations  (44)  and 
(45),  to  be  S(i,  s')\(i,  s').  The  entire  squared  quantity  under  the  path  integral 
can  then  be  identified,  from  equation  (36),  to  be  0(s')k(s')/2ir.  Thus,  0e(s) 


N 

h 

i=  1 


7iF2(s,)g1(i)g2"f^ 


ds' 


k 
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can  be  written  as  a path  average  of  the  band  model  parameter  0 weighted  by 
the  absorption  parameter  k 


0 (s)  = J C c(s')p(s')k(s')0(s')  'Is'  . (46) 

e u(s)ke(s)*/o 


Equations  (40),  (43),  and  (46)  constitute  the  traditional  formulation  of  the 
CG  approximation  for  a statistical  band  model. 

B.  ALTERNATIVE  FORMULATION 

We  now  consider  an  alternative  approach  to  the  CG  approximation  for 
a random  array  of  lines  that  more  nearly  parallels  the  formulation  of  the 
band  model  for  a uniform  optical  path.  The  starting  point  is  the  use  of  the 
CG  approximation  for  an  isolated  line  [equations  (12)  and  (9)] 


W(s)  = 2 Try 


e 


u(s)S  (s)" 

(s)L[2nve(S)_ 


(47) 


As  in  the  formulation  for  a uniform  path,  we  average  this  expression  over  a 
line  strength  distribution  function  P(S).  Since  the  path  is  assumed  to  be 
nonisothermal,  no  explicit  line  strength  appears  in  equation  (47);  only  the 
path-averaged  line  strength  S^(s)  appears.  It  becomes  necessary,  then,  to 
consider  a strength  distribution  function  for  Se>  An  approximate  distribution 
P(Sp)  can  be  established  by  invoking  the  separability  approximation 


By  use  of  the  Schwartz  inequality  for  integrals,  it  can  be  established  that  if 
the  separability  approximations  of  equations  (44)  and  (45)  are  not  invoked, 
then  the  optimum  choice  for  0 (s)  is  some  value  less  than  or  equal  to  the 
path  average  given  by  equation  (46). 
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[equation  (44)]  on  line  strength.  The  path  average  for  the  strength  of  the  i 
line  in  Av  is  [equation  (6)] 


S^U,  s)  = J“  c(s')p(s')S(i,  s')  ds' 

The  use  of  equation  (26)  for  S(i,  s')  yields 

c(s/)p(s/)SM(s/)  e"9(l)/T(s,)  ds'  . (48) 

The  separability  approximation  is  now  introduced  by  assuming  that  T(s')  in 
the  exponential  argument  of  equation  (48)  can  be  replaced  by  a quantity  Te(s) 
that  is  independent  of  s'.  This  replacement  corresponds  to  the  use  of 

-8(i)/T  (s) 

CT^i)  = e 


S (i,  s)  = — : — r 
e ’ u(s) 


L 


and 


-2(s/)  = Sm(s') 


in  equation  (44).  Then, 


-0(i)/T  (s) 

Se(i's)  = SeM(s)  e 


where 


SeM(s) 


c(s')p(s')SM(s') 


ds ' 


(49) 


(50) 
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r 


Regardless  of  the  choice  of  Te(s),  it  is  clear  from  equations  ( 6>)  and  (49) 
that  S (i,  s)  displays  the  same  functional  dependence  with  T (s)  that  S(i,  s) 
does  with  T.  Consequently,  because  the  energy  level  distribution  function 
P(3)  is  independent  of  temperature,  it  is  readily  established  that,  within  the 
separability  approximation,  Sg  is  distributed  according  to  the  inverse 
distribution 


fTe(s) 

s e 

e m 


P(Se)  = < 


0 < S < S ..(s) 

e eM 


S > S . ,(s) 
e eM 


The  application  of  this  distribution  to  W(s)  given  by  equation  (47)  completely 
parallels  the  formulation  for  a uniform  path  given  in  Section  III  and  yields 
the  result 


««(.»]  , 

(51) 

;(s)  = <4£M£>  _ 

0(s) 


wnere 


k(  s) 


Te(l5) 

S (s)-r^ 

e M 0 0 

m 


(52) 


4Te(s) 

3(S)  = 2TT^e(s)  

m 


(53) 
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The  parameter  \g(s)  in  equation  (53)  for  3(s)  is  given  by 


Vs)  g c(s')p(s')SM(s'h(s')  ds'  (54) 

and  is  obtained  from  equations  (6),  (7),  and  (26)  by  again  replacing  T(s  ) with 
T (s)  and  assuming  \(i,  s')  to  be  independent  of  i and  equal  to  the  band  model 
parameter  y(s' ). 

Since  equations  (40)  and  (51)  are  identical  in  form,  it  is  clear  that  the 
assumption  of  equation  (40)  as  the  cx.  ve-of-growth  function  in  Godson's 
formulation  is  equivalent  to  the  use  of  the  separability  approximation  for  line 
strength  and  the  use  of  the  CG  approximation  for  an  isolated  line. 

The  significant  departure  of  this  alternative  formulation  of  the  CG 
approximation  from  Godson's  formulation  is  the  manner  in  which  the  effective 
parameters  k(s)  and  0(s)  are  determined.  In  Godson's  formulation,  the 
effective  parameters  were  introduced  a priori  as  unknowns  to  be  determined. 
The  present  formulation  has  established  a direct  link  between  the  effective 
band  model  parameters  k and  *3  and  the  band  model  parameters  for  isothermal 
paths.  More  explicitly,  equations  (52),  (53),  and  (54)  together  with  equa- 
tion (50)  for  S . ,(s)  give  k(s)  and  i3(s)  in  terms  of  the  modeling  parameters 
em 

S^(s)  and  $3  . These  later  parameters  are,  in  turn,  given  in  terms  of  the 

isothermal  parameters  k(s)  and  6(s)  by  equations  (38)  and  (39).  Caj  rying 
through  the  formulation  with  these  relationships,  we  find  that  k(s),  0(s),  and 
y vs)  can  be  expressed  as 


k(s)  = ke(s) 


rre(s)  Se(sf 


1(s)  t(s) 


B(s)  = 


\c(s)  = 


2-\e(s) 

[T  (s) 
e 

b ( s) 

T(s) 

e 

r c!s- 

^0 

)p(s'  )1 

6(  s)  J 


(55) 
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where  (h)  is  a path  averape  of  the  line  spacing  hand  model  parameter  given 
e 

b> 


1 

u(s)ke(s) 


c(sf )p(s' )k(s' )S(s' ) ds' 


(56) 


and  k (s)  is  given  bv  equation  (43). 

i' 

3 he  weak  limit  form  for  W(s)/‘  with  this  result  for  k(s)  is 


W(s) 


= u(  s )k(  s ) = u(  s)k  ( s ) 
e 


\U) 

Ti (TT 


£(s)  J 


Inasmuch  as  it  has  been  established  (in  Godson's  formulation)  that  the  correct 
weak  limit  form  for  W(s )/L  is  u(s)k^(s),  however,  the  quantity  in  brackets 
must  be  unity.  Thus,  a definition  for  T^Js)  is  established  as 


T (s)  = T(s)f^- 


(57) 


Although  T (s)  does  not  enter  into  the  final  result  of  this  approach  to  the  CG 
e 

approximation,  the  definition  of  T^(s)  by  equation  (57)  plays  ?,n  important  role 
in  the  band  model  formulations  presented  in  Section  V. 

The  final  result  for  this  alternative  formulation  of  the  CG  approximation 
is  the  same  as  Godson's  formulation,  except  that  Pe(s)  is  not  defined  as  a path 
average  of  3(s)  according  to  equation  (46),  but  rather  by 


0 (S) 
e 


2tt\c(  s) 


(58) 


where  \ (s)  is  a path-averaged  line  width  given  by  equation  (55),  and  6 (s)  is 
a path-averaged  line  spacing  parameter  given  by  equation  (56). 
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c. 


DERIVATIVE  FUNCTIONS 


The  computation  of  band  radiance  by  use  of  equation  (19)  requires  the 
transmittance  derivative  d~(s)/ds.  Differentiation  of  equation  (21)  yields 


dT(s) 

ds 


-T(s) 


1 dW(  s) 
6 ds 


(5Q) 


The  equivalent  width  derivative  obtained  by  differentiating  Godson's  result 
for  the  CG  approximation  is 

| d^S-  = c(s)p(s)k(s)y[4xp(s),  r(s)]  , (60) 

where  the  derivative  function  y(x,  r)  is 

y(x,  r)  = (2  - r)^i  + (r  - 1)^  (61) 


with 


djj*l  = iii*i  (62) 

dx  x 


The  nonuniformity  index  r(s)  is 


r(  s) 


gjs) 

FIJI 

e 


(63) 


where  $e(s)  is  given  by  equation  (46).  This  result  is  similar  in  form  to  that 
obtained  for  an  isolated  line  [equations  (13)  and  (14)].  In  particular,  this  re- 
sult can  display  the  same  spurious  behavior  of  y > 1 for  intermediate  x if  the 
nonuniformity  index  r(s)  is  very  large. 
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By  differentiating  the  result  for  W(s)/6  obtained  in  the  alternative 
formulation  for  the  CG  approximation,  we  obtain 

1 ri  vV  ( s ) 

ds — = c(s)p(s)V(s)y(4xt;(s),  p(s),q(s)]  , (64) 


who  re 


y(x„  p,  q) 


q / x 


(65) 


Here,  the  arguments  p(s)  and  q(  s ) are 


P ( s ) ^ ^1-  (66) 

Ve<®> 

and 

6 (s) 

q(s)  r _£ . (67) 

Ms) 

In  the  alternative  formulation,  the  index  p(s)  measures  primarily  the  local 
degree  of  inhomogeneitv  at  path  position  s,  and  q(s)  measures  primarily  the 
local  degree  of  noni  sothcrmality.  The  parameter  r(s)  in  Godson's  formula- 
tion measures  the  combined  effects  of  inhomogeneity  and  nonisothe  rmality. 
l or  an  isothermal  path,  q = 1,  r = p,  and  equations  (61)  and  (65)  become 
identical.  Thus,  the  two  approaches  treat  inhomogeneities  along  the  path  in 
the  same  manner.  The  significant  difference  between  the  two  approaches  is 
the  manner  in  which  nonisothe  rmalitics  are  tieated.  Consider  a homogeneous 
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path  for  strong  abssrption.  Then,  p is  1 and  r - q.  For  large  x,  equations 
(61)  and  (65)  yield,  respectively. 


and 


y(4x,  q) 

(2irx)1/Z 


y(4x,l,q)  ^ 


(2ttx) 


Godson's  approach  yields  a linear  increase  of  y with  q for  strong  absorption. 
The  alternative  approach  displays  a much  weaker  increase  with  q through  the 
func*;on  (2q  - l)/q. 

Tikfc  results  for  either  formulation  of  the  CG  approximation  with  the 
exponential- tailed  distribution  are  identical  to  all  of  the  results  in  this  sec- 
tion, except  that  all  occurrences  of  I(x)  are  replaced  by  1 (x),  and  the  fac*o: 

e 1/2 

4 Chat  appears  with  x^is)  is  replaced  by  it.  In  addition,  dle(x)/dx  = 1/(1  + 2x) 
replaces  equation  (62). 
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V.  DERIVATIVE  APPROXIMATIONS 


A.  INTRODUCTION 

The  approach  to  the  band  models  formulated  in  this  section  is 
fundamentally  similar  to  the  approach  of  the  LS  approximation  for  an  isolated 
line.  The  radiance  equation  [equation  (19)1  for  band  radiation  is  formulated 
in  terms  of  the  transmittance  derivative  dT(s)/ds  rather  than  in  terms  of  the 
transmittance  t(s)  directly.  The  CG  approximations  presented  in  Section  IV 
are  essentially  approximations  to  t(s)  in  that  they  deal  directly  with  the 
quantity  "W(  s)/  6 s -ln~(s).  Expr  e ssions  for  dT(s)/ds  are  obtained  by  differen- 
tiating the  approximations.  In  this  section,  direct  approximations  to  dT(s)/ds 
are  treated.  More  accurately,  because  [equation  (59)] 

d~(s)  _ - v 1 .dWjsj 

ds  " ' 6 ds 


approximations  to  the  mean  equivalent  width  derivative  are  considered.  The 

approach  is  similar  to  the  alternative  formulation  of  the  CG  approximation  in 

that  we  start  with  an  approximate  solution  for  an  isolated  line,  average  over 

a line  strength  distribution,  and  express  the  distribution  parameters  S^  and 

*6  in  terms  of  band  model  parameters  for  isothermal  paths  by  the  direct 
m 

relations  of  equations  (38)  and  (39).  The  results  are  thus  expressed  in  terms 
of  the  nonuniformity  parameters  p(s)  and  q(s).  However,  intuitive  models 
based  on  the  parameter  r(s)  are  also  considered. 


B.  DERIVATIVE  APPROXIMATION 

Where  the  alternative  formulation  of  the  CG  approximation  was  based 
on  the  CG  approximation  to  W(s)  for  an  isolated  line,  we  begin  here  with  the 
LS  approximation  to  dW(  s)/ds  [equations  (9),  (15),  (17),  and  (18)] 


dW(  s) 
ds 


c(s)p(s)S(s) 


2u(s)S  (s)  . 

e i 

dz 

2^eU)  1 + zV(s) 

1 + z2 

(68) 
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Preceding  page  blank 


In  the  CG  formulation,  the  only  line  strength  parameter  in  the  equations  was 
the  path-averaged  strength  Se>  An  average  over  line  strength  was  performed 
by  deriving  an  approximate  distribution  for  Se>  In  the  present  case,  botli  the 
path-averaged  and  local  line  strength  S(s)  appear  in  the  formulation.  In  order 
to  perform  an  average  over  line  strength,  it  is  necessary  to  express  as  a 
function  of  S,  or  vice  versa.  This  is  accomplished  by  introducing  the  line 
strength  separability  approximation,  but  only  for  S(s'  ) as  it  appears  in  the 
path  integral  for  S^(s).  That  is,  equation  (49)  is  used  to  express  S^i,  s)  as 

-G(i)/T  (s) 

Se(i,  s)  = SeM(s)  e . (69) 

The  local  line  strength  parameter  S(i,  s)  is  not  separated  and  is  retained  as 
given  by  equation  (26), 


S(i,  s)  = SM(S)  e 


-6(i)/T(s) 


(70) 


Elimination  of  0(i)  between  equations  (69)  and  (70)  gives  the  desired  approxi- 
mate relationship  between  3 and  S 


•Vs)  = ScM(s)'~ 


T(s)/T  ' s) 
e 


From  equations  (57)  and  (67),  the  temperature  ratio  can  be  expressed  as 


T(s) 


6 (s) 


Ti 


(s) 


q(  s) 
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r 


With  the  definition 


S (s)  can  be  written 
e 


Se(s)  - SeM(s)T| 


Then,  with  equations  (71)  and  (72),  equation  (68)  becomes 


c(s)p(s)SM(s)Tl^/  exp- 


da 


1 + 2 P (S)  1 + Z 


2 


where 


u(s)SeM(s) 

XM  S = 2tty  ( s) 
e 


This  expression  can  now  be  averaged  over  the  line  strength  dis.ribution  function 
for  the  local  strength  S(s).  In  terms  of  the  paramet  Q,  the  distlibution 
( equation  (27)]  is 


P(T|)  = 


0 < T]  < 1 


T1  > i 


The  relation  for  the  mean  equivalent  width  is  then 

c(s,P(s)  sm(s)T^ 


m 


7'  i 

2xM(s)Tlq(s) 

/ exp 

2,  , 2 

dTl 

_ 0 

1 + p ( s)z 

1 + z 


The  preintegral  factor  in  brackets  can  be  identified  from  equation  (29)  as  the 
local  value  of  the  absorption  parameter  M*).  With  the  relationships  and  pro- 
cedures used  in  the  alternative  CG  formulation,  it  can  be  established  that 
XM(S)  = 4xe(s),  where  *e(s)  is  defined  by  equation  (40),  and  3e(s)  is  given  by 
equation  (58).  The  final  result  for  the  derivative  is  then 


1 = c(s)p(s)Ms)y[4xe(s),  p(s),  q(s)j 


(73) 


where 


y(x,  p,  q) 


rw 

' r 1 

/ 

/ exp 

2xTlq 
. 2 2 

dll 

0 

l/o 

i + p z 

m 

dz 


1 + z 


(74) 


The  form  of  equation  (73)  is  the  same  as  those  obtained  in  the  CG  approxi- 
mations, but  the  derivative  function  y(x,  p,  q)  is  different.  In  particular, 
where  the  derivative  functions  in  the  CG  approximations  can  assume  values 


The  application  of  P(Tj)  to  dW(s)/ds  gives  the  mean  value  of  the  derivative. 
Clearly,  the  operations  of  differentiation  with  respect  to  s and  averaging 
over  line  strengths  (T))  are  independent  operations.  Thus,  the  mean  value  of 
the  derivative  is  equal  to  the  derivative  of  the  mean  value  of  W,  i.  e.  , 
dW(s)/df.  = dW(s)/ds. 
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greater  than  unity  for  some  values  of  the  argument,  that  of  equation  (74)  is 
well-behaved  for  all  values  of  the  arguments.  A comparison  of  the  two  CG 
expressions  for  y [equations  (61)  and  (65)]  and  an  approximation  (Appendix) 
to  the  present  function  is  shown  in  Fig.  3 for  the  case  p = 2,  q = 3,  and  r = 6. 

Some  analysis  of  y(x,  p,  q)  as  given  by  equation  (74)  has  been  made  by 


(14) 


Lindquist  et  al.,11’'  and  an  approximate  method  of  computing  y(x,p,q)  is  given 
here  in  the  Appendix."  An  interesting  behavior  of  this  function  not 

predicted  by  the  CG  functions  is  the  asymptotic  form  y~  CONST/x^  for 
large  x when  q > 2.  For  q < 2,  the  expected  relation  y ~ CONST/x 
obtained.  This  feature  is  evident  in  Fig.  4,  where  y(x,p,q)  is  plotted  as  a 
function  of  x for  several  values  of  q and  the  fixed  value  p = 1.  These  curves 
were  generated  with  the  approximate  procedure  outline  in  the  Appendix. 

C r mnoiIIST- SIMMONS  APPROXIMATION 

Lindquist  and  Simmons"11  extended  their  approximation  for  an  isolated 
line  to  the  case  of  a random  array  of  L .rentr.  lines  with  equal  strengths. 
Young"2,  131  further  extended  the  approximation  to  the  case  where  the  line 
strengths  are  distributed  according  to  an  exponential  or  exponential-tailed 
inverse  distribution.  These  latter  extensions  were  made  with  the  assumption 
that  S and  S are  directly  related,  and  that  the  constant  of  proportionality  is 
independent^!  the  line  index  i.  This  assumption  is  clearly  equivalent  to  the 
assumption  that  q(s)  = 1 in  equation  (72).  The  result  obtained  for  the 
exponential-tailed  inverse  distribution  was  given  as 


c(s)p(s)k(s)y[Trx  (s),  r( s )] 
ds  e 


(75) 


In  this  approximate  procedure,  the  exponential-tailed  inverse  line  strength 
distribution  is  used.  If  y(x,P,q)  is  computed  according  to  this  approximation, 
the  factor  4 in  equations  (73)  and  (78)  should  be  replaced  by  w. 
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Intcrcomparison  of  Derivative  Functions  for  Band 
Models.  — CG  approximation  (Godson's  method); 

CG  app roximatior  (alternative  formulation); 

formal  derivative  approximation  [equation  (73)]; 

te,  — - — intuitive  derivative  approximation 
(equation  (78)]. 


where  y(xf  r)  is  the  relatively  simple  expression 


y(x,  r) 


2r(l  -f  x)  + (1  + r‘  ,ly/l  + 2x 
sj  1 + 2x  (r  + 1 + 2x)^ 


(7b) 


The  arguments  xe(s)  and  r(s)  were  given  in  terms  of  3e(s)  computed 
as  a path  average  of  0 [equation  (4 6) | as  is  used  in  Godson's  formulation  of 
the  CG  approximation.  This  choice  for  9e(s)  was  made  intuitively  so  that  the 
strong  absorption  limit  (x  ®)  of  y(x,  r)  given  by  equation  (76)  and  y(x,  r)  given 
by  the  CG  result  of  equation  (61)  would  agree.  If  the  formulation  for  the 
exponential -tailed  distribution  is  again  carried  out  for  q(s)  s 1,  but  without 
the  intuitive  use  of  r(s),  and  with  the  direct  relations  of  equations  (38)  and 
(3b),  the  result  is 


T^lTs^  = c(  s >P(  s >k(  s )y[  7rxe(  P(s)]  > (77) 

where  y(x,p)  is  still  given  by  equation  (’’6),  but  where  t luniformity 

argument  p(s)  is  now  given  by  equation  (66).  Thr  result  equation  (75)  is 
referred  to  herein  as  the  intuitive  LS  approxima  and  that  of  equation  (77) 
as  the  formal  LS  approximation. 

The  intuitive  use  of  r(s)  instead  of  p(s)  in  order  to  match  the  strong 
limit  forms  of  the  LS  and  CG  approximations  raises  the  possibility  of  intui- 
tively replacing  p(s)  by  r(s)  in  the  derivative  approximation  of  equation  (73) 
also,  if  the  result  of  equation  (73)  is  referred  to  as  the  formal  derivative 
approximation,  then  the  form 

T = c(s)p(s)k(s)y[4xe(s),  r(s),  q(s)]  (78) 

with  y( x,  r,q)  given  by  equation  (74)  is  referred  to  as  the  intuitive  derivative 
approximation. 
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D. 


MF  AN-LINE)  APPROXIMATION 


The  following  approximation  is  introduced  to  permit  th°  use  of  t .e 
simple  form  of  the  derivative  function  given  by  equation  (76),  but  to  lessen 
the  d ;gree  of  approximation  inherent  in  the  assumption  q(s)  = 1.  If  the  line 
strength  as  given  by  equation  (26)  is  averaged  over  the  energy  level  distribu- 
tion P(0 ) of  equation  (27),  the  result  is 


S(s)  = SM(s> 


T(s) 

0 

m 


Consequently,  the  relationship  between  the  mean  line  strengths  at  arbitrary 
path  positions  s and  s'  is 


S(  s'  ) 


Vs'  > 


T(s'  ) 

“t nr 


S(s) 


(70) 


The  mean-line  approximation  assumes  that  this  functional  relationship  holds 
for  each  line  in  Av,  thus, 


S(i,  s'  ) 


SM(s'’  T(  s'  ) 

SM(>>  T(sT 


S(i,  s) 


(BO) 


Substitution  of  equation  (80)  into  equation  (6)  for  Se(i,  s)  and  use  of  equations 
(38),  (39),  and  (56)  yield 


S (i,  s) 
e 


(81) 


where  11  and  q(s)  are  defined  by  equations  (71)  and  (6^)f  respectively,  ihis 
approximation  retains  the  linearity  of  Se(i,  s)  with  S(i,  s),  but  defines  the 
constant  of  proportionality  explicitly  in  terms  of  path  properties.  In  particula 


an  inverse  dependence  of  Se(i,  s)  on  q(s)  is  predicted.  Comparison  of  equation 
(»1)  with  equation  (72)  shows  that  the  former  is  obtained  by  factoring  T)^  as 
T|7]^  * and  replacing  ^ by  its  mean  h.e  in  the  interval  0 S ?,  £ 1 


The  elimination  of  S ^j(s)  from  equation  (81)  with  equ„.*ions  (50)  and  (38)  re- 
sults in  the  final  form  for  the  formal  mean-line  approximation  as 


1 dW(  s) 
l ds 


c(s)p(s)k(s)y 


f*rx  (s) 
e 

q(s) 


, p(s) 


(82) 


where  y(x,p)  isi  given  by  equation  (76).  As  for  the  previous  two  approxima- 
tion, we  can  consider  an  intuitive  mean-line  approximation  given  by 


1 dW(  s)  _ 
6 ds 


c(r,)p(s)k(s)y 


ttx  (s) 
q(sj  ' 


r(s) 


(83) 
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VI. 


DISCUSSION  AND  SUMMARY  OF  BAND  MODELS 


A.  SUMMARY  OF  FORMULATION  PRINCIPLES 

The  formulation  of  band  models  for  a random  array  of  lines  requires 
integrations  over  the  spatial  variable  s',  the  spectral  variable  v,  and  the  line 
strength  variable  S.  In  a completely  general  formulation,  these  three  vari- 
ables are  intimately  interrelated  in  a complex  manner.  The  heart  of  nearly 
all  band  model  formulations  for  general  optical  paths  is  the  introduction  of 
approximations  that  separate  the  affects  of  the  variables  such  that  the  orders 
of  integrations  can  be  interchanged.  In  the  models  discussed  here,  the  vari- 
ables s'  and  v are  separated  with  either  the  CG  or  LS  approximations  for  a 
single  line.  The  variables  s ' and  S are  effectively  separated  by  using  the 
line  strength  separability  approximation. 

The  fundamental  feature  of  the  CG  approximations  is  that  they  are 
approximations  to  the  mean  equivalent  width  function  W(s)/6.  The  various 
derivative  approximations,  on  the  other  hand,  attempt  to  make  direct  approx- 
imations to  the  derivative  of  the  mean  equivalent  width  function,  which  is  the 
more  fundamental  of  the  two  from  the  standpoint  of  calculating  band  radiance. 

The  CG  approximation  was  formulated  in  two  ways  ; (1)  by  defining  the 
effective  band  model  parameters  k (s)  and  5 (s)  acco  ruing  to  Godson  s pro- 
cedure  of  fitting  the  model  to  exact  weak  and  strong  absorption  limits  (within 
the  confines  of  the  separability  approximation),  and  (2)  by  defining  these 
effective  band  mode)  parameters  directly  in  terms  of  isothermal  band  model 
parameters  through  relations  that  are  derived  from  the  properties  o)  the 
assumed  inverse  line  strength  distribution  function.  The  two  approaches 
appear  to  be  equivalent  because  they  are  based  on  the  same  assumptions. 
However,  the  results  are  somewhat  different.  The  latter  formulation  d< 
pends  more  on  the  validity  of  the  assumed  line  strength  distribution  model 
than  does  the  former. 
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The  derivative  approximations  were  all  derived  in  a manner  similar 
to  the  second  approach  to  the  CG  approximation.  These  formulations  are, 
the  1,  functions  of  the  inhomogeneity  index  p(s)  and  the  nonisothermality 
index  q(s ).  In  order  to  match  the  strong  absorption  limit  of  one  of  these 
derivative  band  models  (LS)  to  the  corresponding  limit  of  the  first  approach 
to  the  CG  approximation,  an  intuitive  replacement  of  p(s)  with  r(s)  was  intro- 
duced. This  replacement  was  extended  to  the  remaining  two  derivative 
approaches . 

In  all  of  the  models,  temperature  variations  along  an  optical  path  are 
implicitly  accounted  for  in  the  temperature  variation  of  k(s)  and  Ms).  Ex- 
plicit ac'-ount.  is  made  in  all  models,  except  the  formal  LS  model,  by  the  use 
of  either  the  nonisothermality  index  r(s)  or  q(s).  The  intuitive  models 
account  most  for  local  nonisothermalities  by  using  both  r(s)  and  q(s).  Within 
the  heirarchy  of  models,  the  degree  of  variation  with  thermal  gradient  in- 
creases in  the  order  CG,  LS,  mean  line,  derivative  approximation  for  the 
intuitive  models,  and  LS,  CG,  mean  line,  derivative  approximation  for  the 
formal  models. 

The  results  of  the  companion  paper  indicate  that,  for  H^O,  the  intui- 
tive derivative  approximation  yields  the  best  comparison  between  predicted 
and  available  experimentally  measured  radiance  spectra. 

B.  SUMMARY  OF  MODELS 

The  eight  models  considered  here  are  summarized  below  by  extracting 
the  relevant  equations  from  the  text  and  presenting  them  in  groups  that  cor- 
respond to  the  inverse  order  in  which  they  would  be  used  in  making  an  actual 
calculation.  In  all  cases,  the  thermodynamic  variables  c(s),  p(s),  and  T(si 
are  presumed  known  along  tne  general  line  of  sight  that  extends  from  s - 0 
to  s S.  The  bard  model  parameters  k(s),  Ms),  \(s),  and  B(s)  are  pre- 
sumed known  as  a Junction  of  c,  p,  T,  and  thus  s.  All  ol  the  equations  are 
written  for  the  exponential -tailed  inverse  line  strength  distribution. 
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RADIANCE  AND  TRANSMITTANCE  EQUATIONS 


-n 


N = 


ds  , 


llll)  W-)  1 d^^S) 

ds  ' [S)  l ds 


t(s  ) = exp 


W (s ! 

"6 


_1_  dW(s ) 
6 ds 


= c(s  ) p(s ) k(s ) y(s ) 


l£ -)  Iel7rxe(s)^ 


fi 
M s 


S 1 dW(s  ')  , , 
ds 


6 ds 


CG  approximations 


De  rivative 
approximations 


i.  DERIVATIVE  FUNCTIONS 
a.  Godson's  CG 


y(s)  = y[irxe(s),  r(s  ) J , 

dl  (r ) I (T,) 

y(r.r)  = (Z  - r)  — jp — t (r  - 1 ) — = — 

ie(Ti)  = jmT  - i . 

dl  (11)  1 

e = 

di  /fTTH" 


-57- 


b. 


Alternative  CG 


y(s ) = y[Trxe(s),  p(s),  q(s)J  , 


yCH,  p.q) 


c.  Formal  LS 


y(s ) = y[TOe(s  )»  p(s)]  , 


yCH*  p ) 


2P(  1 + T)  T (1  + pZ)  x/1  + 2" 

ynrr  <P  + yi  + 


d.  Intuitive  LS 

The  same  as  formal  LS,  except  with  p(s)  replaced  by  r(s). 

e.  Formal  Mean  Line 


y(s) 


'ttx  (s  ) 
e 

q(s) 


p(s)  , 


y(Tj,p)  = same  as  formal  LS. 


f.  Intuitive  Mean  Line 

The  same  as  the  formal  mean  line,  except  with  p(s)  replaced  by  r(c). 

g.  Formal  Derivative 


y ( s ) =y[;rxe(s),  p(s),  q(s  ) j , 


where  y(”l,p,q)  is  from  the  Appendix. 


h.  Intuitive  Derivative 

The  same  as  the  formal  derivative,  except  w h p(s)  replaced  by  r(s). 

i.  PATH  AVERAGES 


u(s ) 


f c(s  ')  p(s  ')  ds  ' , 

Jc\ 


kc|sl  • ills) 


i—  f c(a  ')  p(s  ' )k(  ' ) ds  ' 

s)  Jn 


6 (s)  = 

e 


i (s  ) ke(s  ) 


/s 

c(s/)  p(s  ' 


) k(s  ')  6(s  ')  ds  ' , 


\„(s ) — _ _ 

c u(s)  k (s)  6 (s) 

e e 


/s 

c ( s ' 

K IM  * - 


) p(s  ')  k(s  ')  Ms  ')  y(s  ')  ds  ' 


6 (s ) 

e 


I 

\ 


u(s)  k (s)  •'o 
c 


f c{s  ')  p(s  ')  k(s  ')  9(s  ')  ds  ' 

Jc\ 


2 Try  (s  ) 
1 e 

T (s) 

e 


The  first  option  for  BJs)  is  u?ed  with  Godson’s  formulation  of  the  CG 
approximation  and  in  the  intuitive  derivative  approximations.  The  second 
option  is  used  with  the  alternative  formulation  of  the  CG  approximation  and 
• he  formal  derivative  approximations 


-59- 


u(s)  k (s) 

><s)  stfr  ■ 

p(s)  \(s)/\e(s)  , 
q(s)  j^(s)/ 6 (s  ) , 

r(s  ) = @(s)/£e(s)  [ 3 e (s  ) from  option  1 above] 


-60- 


APPENDIX 


COMPUTATION  OP  y(x,  p,q)  IN  THE  DERIVATIVE 
APPROXIMATION 


Application  of  the  derivative  approximation  band  model  lequires  the 
computation  of  the  function 


y(x, p,  q)  = 


- f\  (' 

" Jo  |/o 


I 2x^q  ) 
exp  . - — _ 

( 1 + p Z ) 


1 + z 


for  arbitrary  positive  values  cf  the  arguments  x,  p,  q.  The  inner  integral 
may  be  expressed  in  terms  of  the  incomplete  gamma  function 


/ e\  C ^ a-i  'u 

y(a,  ?)  = / u e 


2 2 

by  the  transformation  u = 2x‘rP/(l  + p‘“  z ).  Then,  equation  (A  - 1 ) becomes 


y(x,  p.q) 


I J Z2)l/q  \(~  > ^-7)  dz  . (A -2) 

q(2x)q7T0  1+z  \q  1 + p z / 


Computation  of  y(x,  p,  q)  by  numerical  integration  of  equation  (A-2)  is  com- 
plicated by  the  need  to  evaluate  the  incomplete  gamma  function  and  by  the 
overall  complex  behavior  of  the  integrand.  In  order  to  circumvent  these 
difficulties,  the  following  approximate  method,  which  has  been  found  to  be 
both  sufficiently  accurate  and  efficiently  applicable  is  used. 
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i 


First,  we  write  y(x,  p,  q)  in  the  form 


y(x, p, q) 


- y(x,P, l) 


y(xT,q,  p,  1 ) dT,  , 


(A -3) 


which  expresses  y(\,p,q)  in  terms  of  its  value  at  q = 1 . This  identity  ran 
be  established  as  follows.  For  q = 1,  tue  inner  integration  of  equation  (A-l) 
is  easily  performed  and  yields  the  result 


y(x. p,  1 ) = 


/•cd  ? 2 

4 ^ '—(■  7-77) 


dz 


1 -r  z* 


(A -4) 


Heplacement  of  x with  xTjq  and  integration  of  both  sides  over  0 ^ T)  < 1 yields, 


/ y(x~,q,P,l>  dT(  / fl  1 : / ^ |i  . exp/. 

' */0  !/o  2x~q  I l 1 ^ p“  z2/j 


dT, 


dz 


1 + z 


If  the  inner  integral  is  integrated  by  parts  with 


u 1 * exp 


\ 1 - P z / 


and 


.32 

dv C — — d’  , 

2x~q 
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we  obtain  (for  q £ 1) 


The  integral  in  the  first  term  is  recognized  from  equation  (A-4)  to  be 
y(xf  p,  1',  whereas  the  double  integral  in  the  second  term  is  recognized 

from  equation  (A - 1 ) to  be  y(x,  p,  q).  Thus, 


y(xr,q,p,  1)  dTj  = — i—  y(x,p,  1)  - y-^q  V(x'p'q)  ' 

J0  H 

which,  when  inverted  for  y(x,  p,  q),  yields  and  establishes  the  validity  of 
equation  (A-3). 

Next,  we  approximate  y(x,  p,  1 ) by  the  function  derived  for  the 
exponential-tailed  inverse  line  strength  distribution  function.  This  function 

(13) 

has  been  derived  by  Young  and  is 

, 2p(1  + x)  + (1  + PZ>  >/l  + (A-51 

y(x.p.i)  r ; . ,z 

yrrTx  [p  + y/ 1 + 2x] 


The  equation  that  corresponds  to  equation  (A  - 1 ) for  the  exponential -tailed 
distribution  (with  R -•  °°)  is 


y(x, p,  q)  = 


■HV; 


exp  -T)  - 


2x" 


Z Z i 
+ p z ‘ 


dT| 


dz 

1 + T 


Numerical  evaluation  of  this  integral  is  as  difficult  as  the  direct  evaluation 
of  equations  (A  - 1 ) or  (A-2).  Moreover,  an  expression  comparable  to  equa- 
tion (A-3)  cannot  be  written.  Thus,  use  of  the  exponential -tailed  distribution 
is  delayed  until  equation  (A-3)  is  established. 
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Tin-  approximate  method  involves  the  evaluation  of  y(x,  p,  q)  according  to 
equation  (A  - 4 ) with  y(x,  p,  1)  given  by  equation  (A-5).  The  integral  term  of 
equation  (A-  G is  evaluated  numerically  as  follows.  The  transformation 
u xTjq  gives  the  integral  term  as 


'/(x~q , p,  l ) d 


qx 


M q 


( 1 /q  )-  1 


y(u, p , 1 ) du  . 


For  q 1,  the  integrand  u^*  ql  * y(u,  p,  I ) is  a decreasing  function  with  u. 
The  most  rapid  decrease  occurs  for  small  p [by  inspection  of  equation  (A  - 5 ) J 
anti  large  q.  In  the  limit  p - 0 and  q - !r,  the  functional  form  of  the  inte- 
grand is 


u(  1 *•  2u) 


(A-6) 


Thus,  even  in  the  worst  case,  the  integrand  decreases  no  faster  than  u 
Let  1^  be  the  value  of  the  integrand  at  an  arbitrary  value  Uj.  The  value  of 
u,  > Uj  for  which  the  integrand  has  fallen  to  the  value  fl^  (0  < f < 1)  is  found 
from  equation  (A-6)  to  be 


u 


J_ 

A 


(1  - Zuj) 

T~ 


(A  - 7 ) 


I'he  integration  interval  0 £ u £ x is  divided  into  a series  of  sequential  inter- 
vals according  to  the  algorithm  of  equation  iA-7).  The  last  interval  is  trun- 
cated at  x if  u,  > A nonzero  starting  value  of  u,  = e must  be  used.  Thus, 
within  each  interval  (except  0 £ u £ e),  the  integrand  varies  by  less  than  the 
fa i tor  f.  A two-point  Gaussian  quadrature  is  then  used  to  evaluate  the 


integral  over  each  interval.  * The  general  form  of  each  of  these  subintegrals 
is 


2 


uO/q)-l 


y(u, p,  1 ) du  . 


With  the  transformation 


z = A + Bu  , 


B = 2/(u2  - Uj)  , (A-8) 

A = 1 - u^B  , 


the  integral  assumes  the  standard  Gaussian  form 


B 


u(z ) 


<l/q)-l 


y[u(z),  p J dz 


f 


with 


u(z ) 


(A-9) 


Ihe  two-point  Gaussian  quadrature  approximation  to  the  integral  is 


it 

Ul 


u(/. 


,(»/«!»-  1 


yt  «(/-  j ).  t'  I 


u(/.2)^l/C,)'  1 y lu(/.2>,  p J | 


with  z.  = 10.  57735027. 
1 » ^ 


(A  - 1 0 ) 


With  a starting  value  of  t - 0.01,  and  with  t = 0.  1,  only  15  intervals  are 
needed  to  handle  x as  large  as  2 X 106.  This  compares  with  hundreds  (or 
more)  of  terms  required  in  various  exact  series  expansions  of  equation  (A  - 1 ) 
and  (A-2). 
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With  f - 0.  10,  this  procedure  gives  an  accuracy  of  ~3%  for  the  worst  case 
of  variation  of  the  integrand,  i.e.,  u and~0.7%  for  a nominal  variation 

, -1 

of  u 

The  contribution  to  the  integral  from  the  interval  0 < u < e is  treated 
by  a two-term  series  expansion.  For  small  u,  equation  (A  - 5 ) gives 

2u 

y(u,  p,  1 ) ► ! - YJJ  * 


Thus , 


/ 


( 1 /q , - 1 


y(u, p , 1 ) du — ►qe 


.(i/q) 


1 - 


2e 

(i  + p)  (i  + q) 


(A  - 1 1 ) 


A nominal  value  c =0.01  is  used.  If  x<  0.01,  this  two-term  expansion  is 
list'd  directly  to  evaluate  the  integral  over  0 < u < x. 

A calculational  flow  diagram  for  the  procedure  is  shown  in  Fig.  A-l. 
Note  tha^  the  routine  treats  only  the  case  q>  1.  A similar  analysis  could  be 
made  for  q < 1,  but  this  case  is  not  as  intrinsically  important  as  the  q>  1 
case.  Il  q < 1 . we  set  q 1. 
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LABORATORY  Oi  ERATtONS 


The  Laboratory  Operation*  of  The  A*-o,p»c«  Corporation  it  conducting 
experimental  and  theoretical  inve*tl|atiom  neceeeary  for  the  evaluation  and 
application  of  acientific  advance*  to  new  military  concept*  and  tyttemt.  Ver- 
satility and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  lccomplishment  of  taske  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophyslcs  Labon1---  ' aunch  and  reentry  aerodynamics,  heat  trans- 
fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Leboruory:  Atmospheric  resetions  and  atmoe- 
phcric  optics,  chemical  reactions  iii^olluted  atmospheres,  chemical  reactions 
of  excited  species  ii.  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  rsnging.  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics,  quantum  electronics, 
lasers,  and  electro  optics,  communication  sciences,  applied  electronics,  semi- 
conducting. s oarconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging:  atmospheric  pollution:  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials,  metal 
matrix  compnsitcs  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry:  spacecraft  material*  and  electronic  components  :n 
nuclei  r weapons  environment,  application  of  fracture  mechanics  to  stress  cor- 
rosic*.  and  fatigue  - induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory-  Atmospheric  and  ionospheric  phy  ilcl,  ridu< 
tio  from  tKc  atmosphere 7 density  and  composition  of  the  atmosphere,  aurorae 
and  airflow,  magnetospheric  physics,  cosmic  rays,  gen«r«ti»sn  and  propagation 
of  plasma  wsves  in  the  magnetosphere,  solar  physics,  studies  of  solar  magnetic 
fields,  space  astronomy,  «-rey  astronomy,  the  effects  of  nu<  lear  erplostons. 
magnetic  storms,  and  solar  ictivity  on  the  earth's  atmosphere,  ionosphere,  «.>d 
magnetosphere,  the  effects  of  optical,  electromagnetic,  end  particulate  radia- 
tions in  space  on  space  systems. 
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